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LQR CONTROL OF THIN SHELL DYNAMICS: FORMULATION AND NUMERICAL 

IMPLEMENTATION 

R.C.H. DEL ROSARIO* AND R.C. SMITH* 


Abstract. A PDE-based feedback control method for thin cylindrical shells with surface- mounted 
piezoceramic actuators is presented. Donnell-Mushtari equations modified to incorporate both passive and 
active piezoceramic patch contributions are used to model the system dynamics. The well-posedness of this 
model and the associated LQR problem with an unbounded input operator are established through analytic 
semigroup theory. The model is discretized using a Galerkin expansion with basis functions constructed from 
Fourier polynomials tensored with cubic splines, and convergence criteria for the associated approximate 
LQR problem are established. The effectiveness of the method for attenuating the coupled longitudinal, 
circumferential and transverse shell displacements is illustrated through a set of numerical examples. 

Key words. LQR method, PDE thin shell model, Galerkin approximation, piezoceramic actuators 

Subject classification. Applied and Numerical Mathematics 


1. Introduction. Thin shell models are used to characterize structural phenomena ranging from wing 
vibrations to deformations in a duct due to an adjacent flow field. In full generality, shell equations can be 
used to model structures with variable curvature and irregular geometries. Furthermore, the models can be 
modified to incorporate a large variety of actuators and sensors and can be coupled with adjacent acoustic or 
fluid fields to model coupled systems. In all cases, the displacements of various shell components are coupled 
due to the curvature and geometry. 

We focus here on cylindrical shell models due to their prevalence in applications (e.g., noise control 
in a fuselage or flow control in a flexible pipe). For specificity, we consider surface-mounted piezoceramic 
actuators due to their capability for both sensing and actuating. Their applicability is enhanced by the fact 
that their response is relatively linear at low to moderate drive levels. Furthermore, they are lightweight, 
space efficient, and can be manufactured for various geometries at reasonable cost. We note that while 
the control methods discussed here are specific to surface- mounted piezoceramic actuators, the techniques 
can be directly extended to embedded piezo ceramic actuators as well as electrostrictive or magnetostrictive 
actuators which have been restricted to output levels which are approximately linear. 

While experimental work has demonstrated the potential of piezoceramic actuators in cylindrical shell 
applications [13, 25], their full potential is typically not realized due to limitations in hardware, models, 
approximation methods and control laws. One source of difficulty when approximating shell dynamics and 
designing model-based controllers is due to the mechanisms which provide the models with much of their 
flexibility, namely the coupling between shell components. The physical coupling due to geometry and 
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curvature leads to a significant interlacing of frequencies (e.g., as illustrated in [12], it is common in certain 
regimes to have ten or more linearly independent modes within a 10 Hz range). This places stringent 
requirements on the accuracy of approximation methods and the ability of the controller to avoid spillover 
of energy into unmodeled mechanisms. 

One means of approximation in structural applications is through modal or eigenfunction expansions. 
Because analytic expressions for eigenfunctions can be determined only for very restrictive boundary condi- 
tions (e.g., simply supported) for shells with constant material parameters and hence no actuators, sensors or 
coupling components (normal modes do not exist for shells with piecewise constant parameters [6]), approxi- 
mation techniques or assumed modal expansions must be employed. The coupling and consequent interlacing 
of frequencies places stringent requirements on the accuracy of approximation methods or choice of assumed 
modal basis. The effects of coupling must also be considered when designing a feedback control method 
which employs piezoceramic patches. For example, the tight coupling and subsequent frequency interlacing 
can cause a modal controller constructed solely for resolution of uncontrolled dynamics to experience severe 
spillover into uncontrolled modes upon application of the control input. This motivates the development of 
a more general model-based control method for shell applications. 

In this paper, we consider an LQR full state feedback law for cylindrical shells modeled by Donnell- 
Mushtari equations modified to account for actuator dynamics. An approximate feedback law appropriate 
for numerical implementation is developed in the context of a Fourier /spline-based Galerkin method. The 
Donnell-Mushtari equations characterize the coupling and physics in a large number of thin shell applications. 
Furthermore, the analysis and approximation methods developed in this context can be directly extended 
to more accurate models (e.g., Byrne-Fliigge-Lur’ye model) if the application warrants. The consideration 
of LQR performance illustrates properties of the system and model- based control techniques and facilitates 
investigations regarding issues such as patch number and configuration. While full state measurements 
are not available using current instrumentation, and hence the techniques cannot directly be implemented 
in experiments, they provide an important first step in the design of effective compensators based on state 
estimates calculated using a limited number of observations (see [10]). The consideration of the LQR problem 
also provides a step toward the development of model-based controllers for fully coupled structural acoustic 
and fluid/structure systems involving cylindrical shells. As detailed in [12], the Galerkin approximation 
method employed here is flexible with respect to boundary conditions and facilitates consideration of the 
distributional derivatives which arise through the inclusion of patch contributions in the model. It also 
provides the accuracy required to resolve dynamics in the fully coupled model. From a control perspective, 
the use of the general spline basis eliminates the tendency toward spillover into unmodeled components of 
the system. Furthermore, the approximation method satisfies the adjoint convergence and maintenance of 
uniform decay bounds required for convergence of control gains. This enhances the robustness of the method 
in a variety of applications. 

We note that additional theory concerning stabilization and feedback control of thin shells from a 
PDE perspective can be found in [17, 18, 22]. Furthermore, an alternative framework for establishing 
the convergence of control gains for the abstract operator formulation is given in [16, 19]. Both sets of 
references are theoretical in nature, however, and do not provide numerical methods or simulation results. 
The difference and primary contribution of this paper lies in the comprehensive treatment of the problem 
in a framework which includes model development and well-posedness analysis, a rigorous LQR control 
formulation, incorporation of an efficient numerical method, and numerical examples demonstrating the 
scope of the method. 
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Strong and weak forms of the the thin shell model are summarized in Section 2. These models are 
obtained through modification of the classical Donnell- Mushtari equations to incorporate passive (material) 
and active (actuator) contributions due to the patches. The model is also posed in an abstract form based 
upon sesquilinear forms to provide a framework amenable to establishment of model well-posedness and de- 
velopment of an LQR method. Infinite dimensional LQR methods for systems with no exogenous disturbance 
or systems with a periodic disturbance are discussed in Section 3. In both cases, unbounded control inputs 
result due to the piecewise support of the patches. For the case with no exogenous disturbance, criteria 
guaranteeing the convergence of approximate suboptimal gains to optimal gains for the infinite dimensional 
system are established through application of analytic semigroup theory. The Fourier/spline-based Galerkin 
method of [12] is outlined in Section 4, and numerical examples demonstrating the control method arc 
presented in Section 5. These example demonstrates that through the use of the model-based methodol- 
ogy with general Galerkin approximations, significant attenuation in shell vibrations can be obtained using 
piezoccramic actuators. 

2. PDE Model. Throughout this discussion, we consider a cylindrical shell whose radius is large 
in comparison with the thickness (see [15, 21] for details regarding limits on ratios). In accordance with 
common experimental clamping techniques, the edges are taken to be fixed in the sense that longitudinal, 
circumferential and transverse displacements along with transverse slopes are held fixed. Control is pro- 
vided by surface- mounted piezoceramic patches which are configured in pairs with edges aligned with the 
circumferential and longitudinal axes of the shell. 

To specify the geometry for the model, the longitudinal axis of the shell is taken to be aligned with 
the x-axis as depicted in Figure 1. The length, thickness and radius of the shell are respectively denoted 
by h and R while the region occupied by the middle surface is denoted by To. The shell is assumed to 
have mass density p 8 , Young’s modulus F, Poisson ratio z/, Kelvin- Voigt damping coefficient cp and air 
damping coefficient Note that the assumption of Kelvin- Voigt damping is reasonable and typical for 
many shell materials such as aluminum. Finally, the external surface forces are denoted by q x ^Qe^Qn while 
the displacements of the middle surface in the longitudinal, circumferential and transverse directions are 
denoted by u, v and w, respectively. 



X 

Figure 1 . Thin cylindrical shell with surface mounted piezoceramic patches. 
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It is assumed that s pairs of piezoceramic patches having thickness h pe , density p pe , Young’s modulus 
E pe} Poisson ratio v pe and Kelvin- Voigt damping coefficient co pe are mounted to the inner and outer surfaces 
of the shell. Furthermore, it is assumed that the glue bonding layer provides negligible contribution to the 
structural dynamics (the reader is referred to [8, 10] for details concerning the incorporation of differing 
patch characteristics and bonding layers in the ensuing models). 

2.1. Strong Form of the Modeling Equation. The dynamics of the previously described thin 
cylindrical shell with surface-mounted piezoceramic actuators are modeled by the modified Donnell- Mushtari 
equations 


Rph 


d 2 u 

w 


R 


d N x dN dx 

dx d6 


= Rq x -R^W^S pei (x,6) 


( 2 . 1 ) 


d 2 » dNt 


Rph w- 


d0 




i— 1 


dm, 

de 


-S pei (x,0) 


d 2 w dw d 2 M x 

ph ~d^ + Rp, ~di ~ R dx 2 


1 d 2 Mg o d 2 M x g 

R de 2 1 dxde 


+ N g 


= Rq n -J_ 2 

i - 1 


R 


d 2 (M x ) p 

dx 2 


+ 


1 d 2 (Mp) pex - 

R do 2 


As detailed in [8, 10, 21], equations of this type are obtained through force and moment balancing with 
only low order terms retained. Here M xy M$ ) M$ x and M x $ are internal moments, N X) N$ y N$ x and N x e 
denote internal force resultants, and (M x ) pei , ( M $) pei , (N x ) pei , (A^) pe . designate the external resultants (line 
moments and forces) generated by the i ih patch pair. The indicator function 


S pet (x,0) = S X)2 (x)S 1,2(0), 


where 

5i l2 (x) = | 


1 , X < (Zli + X2i)/2 

0 , x = (x u + i 2 i )/2 

-1 , X > (x u + x 2 »)/2 


Si, 2(0) = | 


1 

0 

-1 


, 0 < (Bn + 0 2l )/2 
, O=($ u + 02i)/2 
, 0 > ( Q\i -h 9 2 i)/2 


delineates the sense of the forces generated by the i th pair. The symmetry of the function arises from the 
property that for homogeneous patches having uniform thickness, equal but opposite strains are generated 
about the point {x l ,0 i ) — ((xu -f x 2 i)/ 2, (0 U + 02*)/2). 

The composite density p consists of the shell density p s in regions devoid of patches and a Unear combi- 
nation of p s and the patch density p pe in regions covered by patchs (see [10] for details). Hence p is piecewise 
constant with discontinuities at the patch edges. 

Expressions for the internal force and moment resultants are derived under the assumption that stress 
is proportional to a Unear combination of strain and strain rate. This yields a model which incorporates 
Kelvin- Voigt or strong internal damping. As detailed in [8, 10], the resultants N x , N x e , Nq X) M x , M x ^, Mq x 
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derived under this assumption are 


( 2 . 2 ) 


and 


( 2 . 3 ) 
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The constant a 3 = (h/2 + /i pe ) 3 - h 3 /8 results from integration through the thickness of the patch while the 
characteristic function 


Xpei(x,0) = 


1 , Xii < X < X2i J 01t ^ ^ < 02i 

0 , otherwise 


delineates the region covered by the i th patch pair with edges at xu,X2 », 0H,02i* Expressions for the 
resultants N$ and Mg can be obtained by replacing e x ,ee : K X) k$ in the expressions for N x and M x by 
£ 0 ,£ X i K 0 , K x, respectively. As detailed in [8, 10, 21], the midsurface strains and changes in curvature for the 
Donnell-Mushtari model are 


du 

1 dv w 

dv 1 du 

dx ’ 

£e= Rd 6 + R ’ 

£z0 = di + Rd 6 

d 2 w 

1 d 2 w 

2 d 2 w 

dx 2 

K0 ~ R 2 d 6 2 ’ 

R dxdO 


Note that for the undamped shell which is devoid of patches, substitution of the strain and curvature 
expressions ( 2 . 4 ) into the resultant equations (2.2) and ( 2 . 3 ) yields the classical Donnell-Mushtari expressions 
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N x = 
Ng = 
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(1 - i/ 2 ) 
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(c.g., see [21]). 

The necessity of including the passive (material) patch contributions in the density and moment and 
force resultants for the combined structure is discussed in [10], While the passive patch contributions may be 
negligible in some cases, it is well documented that for many currently employed control configurations, the 
structural dynamics are sufficiently altered by the presence of the patches so as to yield inadequate model 
fits if the patches are neglected. Experimental results demonstrating this effect for circular plates can be 
found in [7, 10] while independent experimental confirmation for rectangular plates is given in [14]. Further 
analysis for thin shells is provided in [27]. The numerical examples in Section 5 of this paper illustrate that 
the control method is equally applicable when passive patch contributions are neglected or included and 
hence is effective when their passive contribution is significant. 

To characterize the external or active patch contributions, it is typical to start with the assumption that 
the strains generated by a patch are proportional to the applied voltage [8]. Since differing voltages can be 
applied to the outer and inner patches in the pair, we will differentiate between the two with Vn(t) and V& (t) 
used to denote the voltages to the outer and inner patches in the i th pair, respectively. The proportionality 
constant relating the generated strain to the input voltage is designated by d 3 i. As detailed in [8], the total 
external moments and forces generated by the patches are 


x Epe d 3i r/a 2 a 3 \ /<*2 \ T/ 1 

( x) pei - l _ Vpe ' h pe Xpe ' LV2 + 37l) V ‘ 1 V2 m) Vi2 \ 


(2.5) 


(N x \ 



d 3 lU2 

i - 

2/lpe 

~E pe 

d 3 i 

JT Xpi 

n pe 

1 - Vpe 


-E, 


(N e ) pei = T-^-dnXpeiSpeAVii + V i2 ] 

1 l/p e 


where a 2 = {h/2 + h pe ) 2 - h?/A and a 3 = (h/2 + h pe ) 3 - h?/8. When substituted into (2.1), the expressions 

(2.5) provide the input from the patches when voltages are applied. We point out that the characteristic 
functions Xpe t restrict the external patch resultants to the region covered by the i th patch pair. 

The fixed-edge boundary conditions 

dw 

(2.6) u = v — w — — — 0 , x — 0, t 

are used to model the end behavior of the shell. These boundary conditions are appropriate for experimental 
setups in which heavy endcaps prevent edge movement. Note that alternative boundary conditions such as 
simply supported or “almost fixed 11 (see [9]) can be employed if edge movement is suspected. 
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We consider two sets of initial conditions. For general systems, the initial conditions are specified as 

u(0,x) — uo(x) , u(0, x) = v 0 (x) > u>(0,x) = wo(x) . 

To model the long-term dynamics of systems driven by a periodic exogenous force g with period r, we also 
consider the periodic conditions 

u(0, x) = u(r, x) , v(0, x) — v(t,x) , w(0, x) = u/(r, x ) . 


2.2. Weak Form of Modeling Equations. The internal and external moment and force resultants 
(2.2), (2.3) and (2.5) are discontinuous due to the presence of the piezoceramic patches. When incorporated 
in the strong form of the shell model (2.1), this leads to problems associated with the differentiation of 
Dirac distributions. To alleviate these difficulties and reduce smoothness requirements on solutions, it is 
advantageous to consider a weak form of the modeling equations which can be derived from Hamilton’s 
principle (energy considerations). 

The state variables for the problem in second-order form are taken to be y = (u,u, w) in the state space 
H = L 2 (To) x L 2 (T o) x L 2 (Fo). For the fixed-edge boundary conditions (2.6), the space of test functions is 
taken to be V = H^{T 0 ) x H^(T 0 ) x H 2 (T 0 ) where 

i?o (To) = {v € H 1 (To) | *7(0, 0) = V (e, 0) = 0} 

flo(To) = {v € tf 2 (r 0 ) 1 77(0,0) = 77* (0,0) = n ( t , e ) = 77 ,(/,0) = 0} . 


For <I> = ( u,v,w ) and V P = (771,772,773), the H and V inner products are taken to be 

<$,$)„= / phufj 1 dj + / phvrj 2 d~f + / phwrj^ 
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i? 2 ( Ke + l ' pe ' Cl ) 90 2 Upe ^ T dxde 


]} 
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where e X) £&, e x e , k x ,kq,t arc defined in ( 2 . 4 ) and dj = RdOdx. The dependence of the inner product on the 
Young’s moduli is explicitly included in the definition to provide a notation for defining analogous damping 
expressions later in this work. 

The weak form of (2.1), as derived in [10] from energy principles, is given by 


L{ m 


d 2 u 

dt 2 


Vi + RN X ^+ N ex ^ - RqxVi - Rj 2 (N x ) pe , d 7 = 0 

i=l ) 


( 2 . 7 ) 


l I Rph~r} 2 + N e ^ + RN xe ^ - Rq e rj 2 - £(JV e ) pei ^ | dj = 0 


L 


d 2 w _ dw_ r _ d^ri 3 

ris + R ^i a + W. ~ rm -^-r 


1 .. d 2 V3 0 ,, d 2 r] 3 

Mq rt-o 2 Mrf- 


< 96> 2 


dxdO 


-^3 + E 


i=l 


a? + fi (Me)pe - a # 2 


^7 = 0 


for all ^ = (771 , 772, 773) ^ A comparison between ( 2 . 7 ) and (2.1) illustrates that in the weak form, deriva- 
tives are transferred from the discontinuous resultants onto suitably smooth test functions. The problem in 
this form is then amenable to analysis and approximation. 

2.3. Abstract Formulation. To provide a natural framework in which to establish model well- 
posedness and infinite dimensional LQR control methods, we consider an abstract formulation of the model 
based upon stiffness and damping sesquilincar forms. To this end, we define cq : V x V — > € , i = 1, 2 by 


( 2 . 8 ) 


<Tl(*,*) = <(£> *W*,*)v 
< 7 2 ($>, S') = ((c c ,c Dp .)$,tf) v + / fiwr} 3 d'y . 

Note that ((c#, co pe )$, ^)y differs from ((E, E pe )$, ty) v only in that Young’s moduli are replaced by Kelvin- 
Voigt damping coefficients. It can be directly verified 1 

(HI) 

m 

(H 3 ) 


m 

(H 5 ) 


It can be directly verified that the stiffness form <j\ satisfies 

|cti($, #)| < ci\^\v\^f\v , for some c\ G 1 R 

(Bounded) 

Reai( 4 >, <£) > C2 1 ^ | v » f° r some C2 > 0 

(F-EUiptic) 

(ri( 9 , 9 ) = <n(V,Q) 

(Symmetric) 

the damping term <72 satisfies 


|<72($, ^)| < C3 1 ^ | v [ ^ f v j for some C3 € 1 R 

(Bounded) 

Re <72(^5^) > C4\&\y , for some C4 > 0 

(F-Elliptic) . 


Remark 1 . The symmetry of <7\ is dependent upon the choice of shell model and ultimately reflects the 
Maxwell-Betti reciprocity theorem. While the Donnell-Mushtari model yields a symmetric sesquilincar form 
cri, other models such as the Timoshenko shell model will not yield a symmetric form. 
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The control input space is taken to be the Hilbert space U = IR 2s and for \I> G V, the input operator 
B G C{U, V *) is defined by 




i 


r E 

F ° 1=1 


dr)i 

dx 


(N x ) pei ^£ + ^m Pei 


dr)2 

ae 


\ d 2 m i tllf x d 2 m 1 j 
{M x ) pei > dy 

where (■, -) K * v denotes the usual duality product. With the definition g = (l/p/i)[g x , q#, g n ], we can write 
the weak form (2.7) in the abstract variational form 


(2-9) <!/(<). $ >v,v + = (Bu{t)+g(t),ty) v . v . 

To pose the problem in a first-order form amenable for control applications, we define the product spaces 
TL = V x H and V = V x V with the norms 

1(01,02)1 n = 101 lv + 102 \"h 
1(01, 02)lv = 101 Iv + I02|v • 

The state is taken to be z(t) = ( y(t ), y(t)) € H while the product space forcing terms are formulated as 


( 2 . 10 ) 


9(t) 


0 

m 


Bu(t) = 


0 

Bu(t) 


The second-order system (2.9) can then be formulated in the first-order form 


( 2 . 11 ) 


(i(t), A) v . v + a(z(t), A) = (Bu(t) + g(t), A) v . v for A € V 
z (0) = z 0 = ( 2 / 0 , y l) 


where a : V x V — ► (D is given by 


O-(0, 0) = - (02, 01 ) V + <Tl (01, 02 ) + <72(02, 02) 

for 0 = (0i,02 ),0 = (0i,02) € V. As proven in [10, page 109], <r is V continuous and for A > 0, cr(-, •) + 
A (*, ’} n is V-elliptic. From the continuity of <r, it follows that one can define an operator A G £(V, V 4 ) by 

A) = (-4T, A)v*,v . 

To obtain a strong form of the first-order system (2.11), consider the system operator 


( 2 . 12 ) 


doim4 = {(^l, <t> 2 ) G 7i\(j)2 G V, Alfa -f j4 2 02 € H} 




with A\ , A 2 G £(V, V*) defined by 


(Ai<f> l , (f>2)y* y — (01 ? 02) j * — 1,2. 

It should be notated that A is the negative of the restriction to dom.4 of A G £(V, V*) so that a(T, A) = 
(— AT, A) n for T G dorrvl, A G V. A strong form of the abstract system model is then given by 


(2.13) 


z(t) = Az(t) + Bu(t) + g(t) in V* = V x V* 
2 ( 0 ) = Zq . 


The rigorous equivalence of solutions is established through the following theorems. 
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Theorem 1. Under Hypotheses (H1)-(H5) on u\ and 02 , A generates an analytic semigroup T(t) on 
and V* . In terms of this semigroup, the representation 

(2.14) z(t) = T(t)z 0 + f T(t - s)[Bu(s) -h^(s)]c?s 

Jo 

defines a mild solution to (2.13) for z 0 G V* and Bu + g G L 2 ((0, T); V*). Furthermore, this semigroup is 
(uniformly) exponentially stable on V, H and V* . 

Theorem 2. Let z sg denote the semigroup solution to (2.13) given by (2.14) and let z var denote the weak 
solution to (2.9). Under hypotheses (Hl)-(H5), it follows that z ag (zo i J r ) = z var (z^^T) for zo € H and 
T = Bu + £ G L 2 ((0, T); V*). 

Following the convention of [28], we will use the same notation for the semigroups defined on V, H and V* 
since each semigroup is an extension or restriction of the others. Note that donul defined in (2.12) is actually 
dom^v4, the domain of A as a generator of T(t) in Tt. As detailed in Lemma 3.6.1 and Theorem 3.6.1 of [28] 
(see also Section IV. 6 of [24] and Chapter 2, Theorem 5.2 of [23]), the property that A generates an analytic 
semigroup on V, 7t and V* results from the continuity and V-ellipticity of cr. The exponential stability of T ( t ) 
on 7i for second-order systems with strong damping is demonstrated in [4] while the exponential stability 
of T(£) on V and V* in this setting is proven in Lemma 3.3 of [5]. Finally, Theorem 2 is a reformulation of 
Theorem 4.14 of [10] and details can be found therein. 


Remark 2. The previous analysis applies to general V*-valued forces with states having initial values 
zq g V*. In many applications, disturbances are generated by rotating or oscillating components in which 
case it is reasonable to assume that g is periodic. The long-term behavior of the system then satisfies the 
abstract Cauchy equation 


(2.15) 


z(t) = Az(t) -F Bu{t) 4- g{t) 
Z{ 0 ) = z(t) ■ 


in V* where r is taken to be commensurate with all frequencies present in the disturbance g. 


3. Approximation Problem. For implementation purposes, we consider Galerkin approximations to 
obtain trajectories which evolve in finite dimensional subspaces C VcM. To guarantee convergence in 
the subsequent control problem, we assume that the approximation method satisfies the standard convergence 
condition 

(H1N) For any z € V, there exists a sequence z N € such that \z — 1 ►Oasj'V— > oc. 

In this framework, the operator A N : which approximates A is defined by restricting a to x V A 

to yield 

(3.1) {-A n T,A) w - <r(T, A) for all T, A 6 V N . 

For each AT, the Co semigroup on which is generated by A N is denoted by T N {t). The control operator 
is approximated by B N G £(U, V N ) given by 

(3.2) { B n u , A) n = <u, B*A) U for all u G U , A G . 
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Finally, wc let V N denote the usual orthogonal projection of 'H onto V' v which by definition satisfies 

(i) V N T e V N for T € H 

(ii) (V N T - T, A) n = 0 for all AeV w . 

This projection can be extended to V N € £(V* , V ,v ) by replacing the H-inner product (T, A) ?( by the duality 
product (T,A) V . v and considering T € V*. 

The approximate problem corresponding to (2.11) can then be formulated as 

(z N (t), A) n + cr(z N (t), A) = ( B N u{t ) + V N g(t),A) H for all AeV N 

z N ( 0)=V N z o . 


with the solution given by 

z N (t) = T N {t)V N z 0 + f T N (t - s) [B w u(«) + V N g{s )] ds . 

Jo 

In strong form, the finite dimensional system has the form 

z N (t) = A N z N {t) + B N u{t) + V N g(t) 
z N { 0) = V N z 0 . 

This system forms the constraint equations used in Section 4.1. A similar finite dimensional system approx- 
imates the dynamics of the periodic system. 

The exponential stability of the semigroup T{t) was established in Theorem 1. An important issue in 
the subsequent control formulation concerns the uniform stability of the semigroup under approximation 
and this is established in the following theorem from [5, 10] (see specifically Lemma 7.13 of [10]). The use 
of the theorem for our shell system is illustrated in Example 1 of the next section. 


Theorem 3. Assume that the injection i : V H is compact. It is also assumed that the damping 
sesquilinear form can be decomposed as <72 = 5(J\ + <72, for some <5 > 0, where the continuous sesquilinear 
form <72 satisfies for some A € IR 

Re < 7 2 (4, 4>) > ~| 4>\ 2 v ~ ml for all ^ € V . 

Finally, suppose that the operator AfM 2 , where A<i € C{V,V*) is defined by 
compact on V. 

If for some u? 6 IR and M > 1, T(t) satisfies 

I T(t)\ c{n) < Me ut , t> 0, 

then for any e > 0 there exists an integer N e such that for N > Nt, 

\T N (t)V N \ C (H) < Me (u,+£)< , t > 0 
for some constant M > 0 independent of N. 


(a 2 <£, t 7 ) 




V* is 


To obtain matrix representations for the finite- dimensional operators A N ,B N and g N , it is necessary 
to specify a basis for the approximating subspace V N . To exploit the tensor nature of the shell domain 


ll 



To and circumferential periodicity, we employ the basis described in [12] which is constructed with Fourier 
components in 0 and cubic splines in x. Specifically, the component bases for the u,v and w displacements 
are assumed to have the form B Uk (9, x) = e lTn9 B Uj (x),B Vk (6>x) = e lTnd B Vj (x), and B Wk (0, x) = e tme B w .(x) 
where B Uj , B Vj and B Wj are cubic B splines modified to satisfy the boundary conditions. The approximating 
subspace then has the form V N = span{£ Ufc }x span{^ fc }x $pan{#u,*} and displacements are approximated 
through the expansions 

Af« 

u N (t,0,x) = y ^u k (t)B Uk (9,x) 

fc — 1 

(3.4) v N {t,e,x) = y \k(t)B Vk {0,x) 

k= 1 

w N (t,0,x) = ^ ~^w k (t)B Wk {6,x ) . 

1 


The restriction of a to V N and construction of the forcing vectors then yields the matrix system 


Kg 

0 

Kg 

0 


0 


1 1 

■ ■ 

= 

o Ke 

-i\ E ” A C D 


1 1 

MO MO 

£ c 
1 1 

+ 

0 

[«(<)] + 

0 

. . 

0 




V 







m k 


tf*(0) 
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where t ) = [ui(t), - • • , ux u (t), V\ (t), • • • , Vtf v {t), w\ (f), ■ * • contains the M = M u + Af v + Af w 

generalized Fourier coefficients. The 5 patch inputs are contained in the vector u(t) = [u \ (£),*••, u s (t)] T . 
The reader is referred to [12] for details concerning the construction of the mass, stiffness and damping 

matrices Af^, Kg, Kg D , the inputs B^ {t) and the initial conditions zg>zg. 

Multiplication by the inverted mass matrix yields the Cauchy equation 


(3.5) 


z N (t) = A N z N (t) + B N u(t) + g N (t) 

Z N ( 0 ) = z» , 


where z N (t) = [^(i), t^(f)] T e with N = 2AT. Note that A N and B N are the matrix representations 
for the operators *4^ and B N in terms of the Fourier/spline basis while g N (t) is the projection of g(t ). 
We point out that the notation z N (t ) designates both the time-dependent generalized Fourier coefficients 
in IR^ and approximate solutions in The specific usage is indicated by the context and the use of 

finite dimensional operators or the corresponding matrix representations. Similarly, the control input u(£) 
should not be confused with the longitudinal displacement u(£,x,0). Equation (3.5) with g N (t) = 0 forms 
the constraint for Example 2 in Section 5. 

4. LQR Control Problem. As noted in (2.13) and (2.15) of the last section, the thin shell model can 
be posed in the abstract forms 


(4.1) 


z(t) = Az(t) + Bu(t) + g(t) 
z( 0) = z 0 
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or 


(4.2) 


z(t) = Az(t) + Bu(t ) + g(t) 
z ( 0 ) = z(t) , 


in V*, depending upon the periodicity assumptions on the disturbance and state. In this section, infinite 
dimensional and approximating finite dimensional LQR control methods are developed for the system (4.1) 
with g = 0 and system (4.2) with periodic disturbance g and periodic end conditions. 


4.1. No Exogenous Input. The infinite horizon optimal control problem for systems with no exoge- 
nous input g can be posed as follows: determine the control u which minimizes the quadratic cost functional 

(4.3) J(u,2o) = jf {\Qz(t)\^ + \n^ 2 u(t)\u}dt 

subject to 


z(t) = Az(t) + Bu(t) 


z(0) - z 0 . 

The nonnegative, self-adjoint operator Q ~ T>*T> e C(H) can be specified to weight state components while 
the positive, self-adjoint operator 71 = (7 Z 1 ^ 2 ) 2 E £{U) serves to weight the control. In this case, 1Z acts as 
a soft constraint to prevent excessive voltages to the patches. 

To guarantee the existence of a unique Riccati solution and an exponentially stable closed loop semigroup, 
it is sufficient to require that (*4, B) is stabilizable and (*4,P) is detectable. Under these conditions, it is 
established in [5, 10] that the optimal control minimizing (4.3) is given by 


(4.4) 


u(t) = —7 Z l B*Ilz(t) 


where II solves the algebraic Riccati equation 

(.4*11 + IL4 - TlBK- l B*Il + Q)z = 0 for all 2 e V 

and z(t) = S(t)zo . Here S(t) is the closed loop semigroup generated by A — B7Z~ 1 J3*U. 

The corresponding finite dimensional problem concerns the determination of suboptimal controls u which 
minimize the functional 

J N (u,z 0 ) = £ {iq'V'WR + 

subject to 

i Ar (f) = ^ z ;v (t) + e N «(<) 

z w (0) - V N z 0 

(see (3.3) with g = 0). The following theorem specifies conditions which guarantee the existence and 
convergence of these suboptimal controls to the control (4.4) for the original infinite dimensional system. 
This is Theorem 7.10 of [10]) and details can be found therein. The application of this theorem to the shells 
system is illustrated in Example 1 which follows. 


Theorem 4. Assume that the injection i : V c — > H is compact. Let the sesquilincar form a associated 
with the first-order system (2.11) be continuous and V-elliptic. Assume that the operators A,B,T> satisfy: 
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(A, B) is stabilizable and (.4, T>) is detectable where B € C{U, V*) is unbounded. Consider an approximation 
method which satisfies (H1N). Finally, suppose that for fixed TVo and N > N 0 , the pair (*4^, B N ) is uniformly 
stabilizable and [A N ,V N ) is uniformly detectable. 

Then for N sufficiently large, there exists a unique nonnegative self-adjoint solution 11^ G £(V*, V) to 
the N th approximate algebraic Riccati equation 

+ u n a n - n N B N n- l B N 'n N + = o 


in V" There also exist constants Ms > 1 and 1^3 > 0 independent of N such that the closed loop semigroup 
S N (t) — -b n ti 1 b n sa ti s fies 


S N (t)\ vN <M z e-^ , t>0. 


Moreover, the convergence of the Riccati and control operators 

U N V N z-^Uz inV for every z € V* 

b n 'u n v n - zrnl ->o, 

I c(n,u) 


as N— >oo, is obtained. 


Example 1. 

We consider in this example a shell with constant parameters i/, cp. Such a case would arise if 
modeling a homogeneous shell or a shell in which the variance of material properties across regions with 
actuators is negligible. The sesquilinear forms for this model are specified in (2.8). Due to the constant 
coefficients, cr 2 can be written as a 2 = 5<t\ + <r 2 where S — ^ and a 2 (<I>, S') = /i / Fo wrjsdy. It follows 
immediately that 



for all (j> E V. The boundedness of the operator A 2 generated by <r 2 follows directly from the boundedness of 
d 2 . Furthermore, it is noted that A^^ 1 € C(V*,V) can be written as an operator on — > V by A^ 1 = 
where the injections i : V H,i* : H V* are compact. Thus A^[ l is compact on V which implies that 
i4j"M 2 is compact on V since it is formed from the product of compact and bounded linear operators. 
Finally, the exponential stability of T(£), the stabiliz ability of (A,B) and the detectability of (A, V) are 
guaranteed by Theorem 1. The hypotheses of Theorem 3 are then satisfied for this system and one obtains 
uniform bounds on the approximating semigroups. The convergence of the Riccati and control operators is 
then obtained from Theorem 4. 


Remark 3. We note that an alternative means of establishing the well-posedness of the closed loop control 
problem and convergence criteria for suboptimal control gains is through the utilization of uniform analyticity 
conditions for the semigroup T(t) generated by A as discussed in [16, 19]. For some systems with unbounded 
input operators, such conditions are more readily verified than the requirements of uniform stabilizability 
and detectability. 
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4.2. Periodic Exogenous Input. The second case under consideration is that in which the exogenous 
force g models periodic disturbances such as noise generated by rotating engine components (e.g., propellers 
or turbines) or periodic electromagnetic sources. Under the assumption that g is periodic in time, the 
functional to be minimized is 

Mu) = l £{\Qz(t)\H + \nV 2 u{t)\v)dt. 

subject to (4.2). The endpoint r is chosen to be commensurate with all frequencies present in g . The 
theory for this case is less complete than that for systems with no exogenous input and is currently limited 
to bounded control inputs B. Under the assumption that (.4, B ) is stabilizable, (.4, T>) is detectable and 
g e L 2 (0,r; W), it is verified in [11] that the Riccati equation 

.4*11 4* IL4 4 UB7Z~ l B*n + Q = 0 

has a unique solution. Furthermore, if r denotes the r-periodic solution of the adjoint or tracking equation 

r(t) = -[A - B'7l~ 1 B*U]*r(t) + Ug(t) 
r( 0) = r(r) 

and z is the closed loop solution of 

z(t) = {A- B7Z l B*U}z(t) - BK~ l B*r(t) 4 g{t) 
z(0) = z(t) , 

then the optimal control is given by 

(4.5) u(t) = B*[Hz(t) — r(t)] . 

For implementation purposes, the suboptimal control is computed through approximation techniques anal- 
ogous to those described in the previous section for the system with no exogenous input. The theoretical 
extension of the periodic theory for unbounded input operators B is currently under investigation. The 
practical efficacy of the method in this regime is demonstrated through numerical examples in the next 
section. 

5. Numerical Examples. To illustrate the previously described control laws, we consider in this 
section several numerical examples which demonstrate the attenuation of shell vibrations using surface- 
mounted piezoceramic actuators. The first example employs the feedback law described in Section 4.1 and 
demonstrates the capability of the PDE- based controller for attenuating transient vibrations. The second 
and third examples demonstrate the capability of the period feedback law developed in Section 4.2. The 
shell is driven by a single frequency source in Example 3 while the attenuation in the case of a multiple 
frequency input is demonstrated in Example 4. 

For all three examples, the exogenous disturbance was constructed to model a periodic noise source 
localized near the shell center at 9 = 0 and 6 = 7r. To accomplish this, the surface forces were taken to be 

= ^q x (t)e- 2O{x -* )2 e(0) 

( 51 ) qe = 0 

q n (t,X,B) = q n (t)e- 20 ^M e{ 0) 
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where 


f 1 - e- {e ~ n/ y /2 , 0 < 0 < w 

0 ( 0 ) = { . 

\ 1 - e-V-W/* , n <0 < 2 -k . 

Here x = £[2 denotes the axial center and q x {t) and q n {t) designate the periodic temporal components to be 
specified in the examples. The spatial components of the Force distributions are illustrated in Figure 2. The 
normal force in the x and 0 components can be interpreted as modeling periodic pressure sources adjacent 
to the points (x,0) = (x,0) and (x,0) — (x,tt). 

Twelve pairs of piezoceramic patches were employed as actuators (hence s = 12 in (2,1) and (2.7)) with 
two sizes considered; the large patches had dimensions of x = 0.2,0 = n/ 3 while the smaller patches had 
sizes of x = 0.1,0 = 7r/6. As indicated in Figure 3, the patch pairs were configured in sets of three along 
the lines 0 = 0, 7r/2, 7r, 37 t/ 2. We reiterate that the voltages to the inner and outer patches in each pair 
were independent so as to provide capabilities for generating both inplane forces and bending moments. The 
dimensions and material properties of the shell and patches arc summarized in Table 1. 

To resolve the uncontrolled and controlled shell dynamics, the Fourier index m was taken from -4 
to 4 for a total of 9 Fourier coefficients. The axial components of the longitudinal and circumferential 
displacements were approximated using 13 modified cubic splines while 11 splines were employed in the 
transverse displacement (two degrees of freedom are lost to accommodate the additional zero-slope boundary 
conditions indicated in (2.6)). This yielded M u = = H7 t Jsf w = 99 in (3.4) and a total of N = 666 

coefficients in the ODE system (3.5). 

The important issues of patch number, placement and size fall outside the scope of this work and are 
under current investigation. The number of patches was fixed in these examples to facilitate comparison and 
demonstrate the capabilities of the control method under uniform conditions. 
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Figure 3. Configuration of large patch pairs 1-6 and 10-12. Pairs 7-9 are centered along 6 = 7r. Observation 
lines L\ = {(x,0)[O < x < i^Q = 7t/6}, L<i — {{x,ff)\x = 3^/4, 0 < 6 < 27r} and observation point 
Pi = (x,0) = (3^/4, 7 t/ 32). Configuration of small patches is similar. 



Dimensions 

Parameters 


h = .00127 m 

p = 2700 kg/m 3 


R = Am 

E = 7.1 x 10 10 N/m 2 

Shell 

t = 1 m 

cp = 1.47 x 10 5 Nms 
v — .33 

p = 58.97 Ns/m 2 


hpe = .0001778 m 

p pe — 7600 kg/m 3 
Epe = 6.3 x 10 10 N/m 2 

Patches 

Centers (x, 6): (.25,0), (.5,0), (.75,0) 

cd p ' = 1.7 x 10 5 Nms 


(■25.Tr/2), (.5, -rr/2), (.75 ,tt/2) 

v pe = .31 


(.25, tt), (-5, 7r), (.75,7 r) 

(.25,3tt/2), (.5,377/2), (.75,3tt/2) 
Dimensions Large Patch: x : 0.2 cm, 0 : 7r/3 
Small Patch: x : 0.1 cm, 0 : n/6 

d 31 = 190 x 10 -12 m/V 


Table 1 . Dimensions and physical parameters for the shell and patches. 


Example 2. 

We illustrate here the feedback law developed in Section 4.1 with no exogenous force. The shell was 
initially driven for t = .01 seconds by the forces (5.1) with 

f lOOOTrt , 0 < t < .01 

Qx{t) = q„(t) = < 

{ 0 , .01 < t < .05 

at which point the forces were terminated and the shell was allowed to freely vibrate. The open loop response 
at the point pi as well as the rms displacement values along the lines L\ and are plotted in Figures 4-6 
(see Figure 3 for the specific locations of pi, L\ and L 2 ). The rms plots illustrate standing waves in all three 
components of the displacement. 
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(a) (b) 



(c) 


Figure 4. Uncontrolled and controlled shell displacements at the point p\ with control initiated at t = .01 
using small patches as specified in Table 1; (a) longitudinal u N {t,x,0), (b) circumferential v N (t,x,0 ), (c) 
transverse w N (t, x,0); (uncontrolled), — — (controlled). 


Under approximation, the feedback law from Section 4.1 has the form 

u(t)^~R~ 1 (B N ) T Il N z N (t) 

where 11^ and z N (t) respectively solve 

(5.2) (; 4 JV ) r n JV + u n a n + u n b n r- 1 (b n ) t u n + q n = o 

and 

z N (t) = [A n - B Jv r , (B ff ) T n w ] z N {t ) 

(5.3) 

Z N (to ) = 20 

(see (3.5) for the component matrices A N and B N ). For this example, we neglect the passive patch contri- 
butions and construct the system matrices using constant parameters p, E, v, c D . As indicated in Example 1, 
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(a) 


(b) 



(c) 


Figure 5. Uncontrolled and controlled shell displacements at the point p\ with control initiated at t — .01 
using large patches as specified in Table 1; (a) longitudinal u N {t^x^O), (b) circumferential v N (tyX,0), (c) 
transverse w N {t,x, 9)\ (uncontrolled), — (controlled). 


the convergence and exponential closed loop stability of control gains for this case are provided by Theorem 4. 
The passive patch contributions are then included in the final two examples. 

The state and control weights were taken to be 


' V(di) 

■ K% 

V(di) _ 



i = 1,2,3 


R = V( Tl ) , i = 1, ■ * - , 24 

where M ^ and denote the mass and stiffness matrices and T>(di),'D(d l ), T>(ri) are diagonal matrices 
whose i th component is di,d l ,r l , respectively. The simulations reported here were computed with the values 
di = di = 10 13 and r* = 50. Large values of di,di are not uncommon for such systems (e.g., see [2]) and 
simulations have been run with values as large as di = di = 10 20 without experiencing degradation due to 
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conditioning. Finally, the choice of a weighted mass matrix for the design matrix Q N is motivated by energy 
considerations as detailed in [3]. 

The closed loop dynamics of the shell were computed by numerically integrating the ODE system (5.3) 
over the time interval [.01, .05]. The time history of the component displacements at the point pi for the 
smaller (.1, n/6) patches are plotted in Figure 4. Corresponding time trajectories with larger (.2, tt/ 3) patches 
are plotted in Figure 5 with corresponding rms values along the lines L\ and L 2 plotted in Figure 6. The 
time histories illustrate that while substantial control is obtained with the smaller patches, the attenuation is 
enhanced significantly through the use of larger patches. In this latter case, the vibrations in all components 
are completely attenuated by time T = .015 seconds. The rms plots demonstrate that strong attenuation is 
attained across the shell including regions not covered by patches. This illustrates both the capabilities of the 
patches as well as the issues associated with patch size. The larger patches will be employed throughout the 
remaining examples. Questions regarding patch size, placement and number are under current investigation. 



circumferential fine L 2 ; (uncontrolled), — (controlled). 


Example 3. 

The remaining examples illustrate the feedback law of Section 4.2 which accommodates a periodic 
exogenous disturbance g. For this example, a single frequency 500 Hz temporal input 

q x (t) = q n {t) = 100sin(10007rt) 
was employed in the exogenous forces (5.1). 

The approximate control law for this case was computed in the following manner. The adjoint solution 
r N (t) was computed by numerically integrating the system 

f N (t) = - - B N R- 1 (B N ) T n N ] T r N (t) +U N g N (t) 

r N (r) = 0 
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from the future time r back to the present (the Riccati solution satisfies (5.2)). The optimal control and 
closed loop solution were then obtained via 

(5.4) u{t) = -Br\B N ) T [n N z N (t) - r N (0] 

and 

z N {t) = [A n - B n R-\B n ) t n w ] z N (t) - B N R‘ 1 (B N ) T r N (t) + g N (t) 
z N { 0) = 0. 

To illustrate the effects of passive (material) patch contributions, we consider two cases. In the first, 
passive contributions are neglected, as in Example 2, and system matrices are constructed using constant 
parameters p , E 1 , v , cd * In the second case, all passive contributions are retained in the resultants (2.2), (2.3) 
and density expressions, and hence are included in the system matrices A N and vectors B N ,g N through the 
inverted mass matrix. 

The uncontrolled and controlled dynamics at the point p\ for the first case are plotted in Figure 7 while 
corresponding plots for the system which includes passive properties are given in Figure 8. In both cases, 





(c) 


Figure 7. Uncontrolled and controlled displacements at the point p\ without passive patch compo- 
nents included in system matrices; (a) longitudinal u N (t,x, ^), (b) circumferential v N (t,x y 0), (c) transverse 
w N (t,x,0); (uncontrolled), (controlled). 
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the time intervals are restricted to [0, .05] to illustrate the system dynamics. A comparison indicates that 
the inclusion of passive or material properties has a significant effect on both open and closed loop dynamics 
due to the additional density, stiffness and damping provided by the patches. Hence for many systems, 
it is necessary to incorporate these components to attain an accurate model and a feasible model-based 
control method (see [10, 14] for experimental evidence demonstrating the necessity of including passive 
patch components). Prom a control perspective, it is noted that for this patch configuration, transverse 
vibrations are reduced significantly throughout the shell with attenuation levels of approximately 80% at 
the point p \ . Similar attenuation is noted in the longitudinal and circumferential components. The spatial 
nature of the uncontrolled and controlled shell in this regime is indicated by the rms values along the axial 
line L\ and circumferential line L 2 as plotted in Figure 9. These plots illustrate that in spite of the persistent 
exogenous force, significant attenuation is attained in all components of the displacement. 



Time (seconds) 

(c) 

Figure 8. Uncontrolled and controlled shell displacements at the point p\ with passive patch compo- 
nents included in system matrices; (a) longitudinal u N (t, £,0), (b) circumferential v N (t,x,6), (c) transverse 
w N (t,x,9)\ (uncontrolled), — (controlled). 
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Figure 9. Root mean square (rms) displacements with passive patch components included in system; (a) 
axial line L\ and (b) circumferential line L 2 ; (uncontrolled), (controlled). 


Example 4. 

This final example demonstrates that the performance of the controller is not degraded by a larger number 
of frequencies in the exogenous disturbance and hence structural response. In this case, we employed 

q x {t) = q n (t) = 50 [sin(3207rt) + sin(4407rt) -f sin(6007rt) + sin(6607rt) + sin(8207rt)] 

so that the exogenous disturbance has five frequencies ranging from 160 Hz to 410 Hz. 

The closed loop solution computed using the feedback law (5.4), with the five frequencies included, is 
compared with the uncontrolled solution in Figure 10 and 11. The applicability of this law on the time 
interval [0, .05], which is one half the fundamental period of the 330 and 410 Hz force components, is due to 
the trigonometric nature of these components. The temporal trajectories at the point pi reflect the multiple 
frequencies in the exogenous disturbance and illustrate that the feedback law (5.4) is highly effective for such 
broadband responses. This observation is reinforced by the rms plots in Figure 11 which indicate significant 
reductions along the lines L\ and L 2 . 

The voltages given by (5.4) for patch pairs 1-6 are plotted in Figure 12 (the outputs from the remaining 
patch pairs arc identical due to symmetry). Due to the axial symmetry of the transverse force q nj the voltages 
to patch pairs 1 and 3, 4 and 6 are nearly identical. It is also noted that the voltages to the central patch 
pairs 2 and 5 are larger than those to surrounding pairs to accommodate the larger axial force distribution 
delivered to the shell center. The predominance of the circumferential force near 6 = 0, 7r leads to larger 
voltages to patch pairs 1-3 than 4-6. A comparison of the voltages to the inner and outer patches in each 
pair first illustrates that the voltages are nearly diametrically out-of-phase which leads to the generation 
of significant bending moments (see (2.5)). As indicated in Figure 13, the voltages also exhibit magnitude 
differences and hence also generate inplane forces in addition to the bending moments. The small magnitude 
of the inplane voltages, as compared with the total voltages plotted in Figure 12, illustrates the degree to 
which bending motion dominates inplane motion for this force configuration. Finally, it should be noted 
that the voltages to the patches are less than 80 V rms which is physically reasonable for patches having 
the specified thickness. The combination of these effects provides the controller with significant vibration 
attenuation capabilities for shell systems excited by general inputs. 
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(a) (b) 

Figure 11. Root mean square (rms) displacements; (a) axial line L\ and (b) circumferential line L 2 ; 

(uncontrolled), — (controlled). 


24 
















Figure 13. Sum of inplane voltages to patch pair 1. 

6. Conclusions. This work summarizes issues concerning the formulation and numerical implemen- 
tation of a PDE-based LQR control method for thin cylindrical shells with surface-mounted piezoceramic 
actuators. From a physical perspective, such actuators are popular since they are lightweight, relatively 
inexpensive and provide broadband control inputs. When mounted in pairs on a cylindrical shell, they can 
be used to generate both inplane forces and bending moments which provides the capability for controlling 
the coupled longitudinal, circumferential and transverse vibrations of the structure. 

Controller design for thin cylindrical shells is significantly more difficult than for plates or beams due 
to the complexity of the structures and the inherent coupling of component displacements due to curvature. 
The models, numerical methods and control laws must incorporate this coupling to be fully effective in 
applications involving shell- like structures. 

In this work, Donnell-Mushtari equations modified to incorporate passive and active patch contributions 
were used to characterize the shell dynamics. This model incorporates the primary coupling mechanisms 
while remaining sufficiently simple to permit initial analysis of the control methods. As detailed in [10, 15, 
21, 26], the model can be directly extended to attain higher accuracy if the application warrants. 

A Galerkin method utilizing bases constructed from tensored Fourier polynomials and modified cubic 
splines was used to approximate the system dynamics. As discussed in [1, 20], an important issue when 
approximating shell dynamics concerns the development of methods which avoid shear or membrane locking. 
One manifestation of locking is the existence of model dynamics which are incorrectly approximated by 
the numerical method. The use of a numerical method which exhibits locking can lead to a loss of control 
authority and potential controller destabilization if the approximations are sufficiently inaccurate. As detailed 
in [12], the Galerkin method employed here avoids locking and provides highly accurate approximate solutions 
with moderate discretization levels. 

The LQR control method discussed here accommodates the unbounded input operators due to the 
discrete nature of the patches, A well-posedness and convergence framework is obtained through the abstract 
formulation of the control problem in terms of sesquilinear forms and associated linear operators. Because 
the LQR method requires full state knowledge for computation of control voltages to the patches, it cannot 
be directly implemented in experiments. It does, however, provide a means for quantifying optimal patch 
capabilities and performance and a setting in which to investigate questions concerning patch placement, 
number, et cetera. It also provides the system and input matrices necessary for later construction of a state 
estimator and compensator which can be experimentally implemented [10]. 


26 



The numerical examples demonstrate the performance of the method for a variety of exogenous dis- 
turbances. Specifically, they illustrate that through the incorporation of coupling mechanisms in the PDE 
model and hence controller, optimal patch voltages are determined which significantly attenuate all three 
components of the shell displacement. Furthermore, by modeling the global shell dynamics and patch inter- 
actions through coupled PDE and constructing the control laws in terms of these PDE, significant reductions 
in both transient and periodic displacements are obtained throughout the shell, including regions which are 
devoid of patches. This illustrates the optimal attenuation possible for the given patch configuration and 
operating conditions. 

An important issue which lies beyond the scope of this paper concerns the optimal number, size and 
placement of patches for a given set of operating conditions. These issues transcend the specified control 
method described here and must be addressed when designing any piezo ceramic-based control method for 
the system. The advantage of the PDE-based method for addressing these problems lies in the incorporation 
of basic physical principles in the models and the utilization of these principles in the control laws. 
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